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non-dynamical scalar field we derive the color dielectric generalization 
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1 Introduction 



There exists several attempts to understand dynamics of gluonic fields in 
non-perturbative region. Among them color dielectric models seem to be 
especially promising pQ. Although none of them describes the full range 
of observable phenomena in a satisfactory way, they are simple enough to 
carry out an extensive study of classical solutions and provide guidelines for 
constructing correct effective theory for non-perturbative QCD. 

A class of color dielectric-like models was proposed by Dick in the late 
90's. It was inspired by low-energy limit of certain string theories. In the 
simplest case, with one scalar field coupled non-minimally to the SU (2) gauge 
field, Dick models reproduce a broad family of quark-ant i- quark confining 
potentials. As it has been pointed out in our previous papers |2] one can 
easily fit potentials originating from these models to phenomenological ones, 
obtained from heavy quarkonia spectroscopy. 

Correct potential between quarks is of course only one in the whole fam- 
ily of factors to be considered and reproduced by an effective gluodynamics 
theory. Acceptable model should also exclude colored states from physical 
spectrum, reproduce energy-momentum tensor trace anomaly as well as pre- 
dict existence of massive particles build entirely of gluonic fields - glueballs. 
The last issue is particularly challenging and was addressed by several au- 
thors in the last few months (cf. eg. [31, 111)- Glueball-like interpretation 
has been also proposed for dilaton field in the Dick model IE]- However, 
so far status of the glueball state in the framework of color dielectric models 
is not clear. 

The most natural generalization of the Dick model emerges when dilaton 
field is replaced by a pair of two scalars. Theory of this kind was recently 
studied by Bazeia at al. [7j in the context of confining potential in a system 
of two planar domain walls. They concluded that a theory with two scalars 
is generally better suited to description of anti-screening effects. As we will 
show below, color dielectric model with two scalars possesses also other strik- 
ing features. Surprisingly observed effects are not a simple superposition of 
well-known results for one scalar. In particular, non-minimal two scalar- 
gauge coupling has glueball as well as meson interpretation and simultane- 
ously gives linear confining potential for external sources. For non-dynamical 
scalar fields it is also equivalent to the well-known Pagels-Tomboulis model 
and its generalized version, which correctly reproduces the trace anomaly in 
the language of effective models. 

The plan of this paper is the following. In the next section we introduce 
and shortly justify the most general form of the model discussed. In section 
El standard analysis of Coulomb problem for external source in the case of 
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massless scalar field is carried out which eventually leads to confining be- 
havior. It also allows us to obtain bounds on the parameters and rule out 
the whole family of models which do not lead to correct confining potential. 
Afterwords we deal with magnetic monopoles and discuss their possible ap- 
plication to glueball physics. In section H] we present similar discussion for 
massive scalars. Finally, non-dynamical scalars are studied in section 



2 The model 

In the present paper we focus on the following density of the Lagrange func- 
tion 

(1) 

where the coupling between scalar and gauge fields is chosen to be 




Here A is a dimensional constant: [A] = cm -1 and a, b are positive parame- 
ters. A particular form of the potential V(<p,ip), fixing asymptotic value of 
the scalar fields will be specified later. 

We would like to notice that all analytical solutions obtained below can 
be also found for the more general coupling function 



<t(0,<0) = 




Due to the fact that for — > oo, ip — > oo the gauge part of the model 
takes the standard Yang-Mills form i.e. a — > 1, the short range behavior 
of the fields changes. For instance, the electric field falls as r -2 for r — ► 0. 
One can check that in the presented solutions the additional, well-known 
Coulomb term will appear. However, such form of the coupling does not 
modify the most interesting long range behavior of fields. Thus we will 
neglect the denominator in the coupling function that is restrict discussion 
to the dielectric function given by (j2J. 



3 The massless scalar fields 

Firstly we analyze the simplest case where the scalar fields are massless. Ad- 
ditionally, the potential term is neglected. Due to that, as we will show it in 
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the next subsection, the asymptotic values of the scalar fields can be arbi- 
trarily large. These values correspond to the so-called dilaton and modulus 
charge [E]. The Lagrange density has the following form 

L = -\ (^) F^F^ + \d,<j>d^ + \d^d^. (4) 

The pertinent equations of motion read 

D„ [(-) I -| r (0) 




and 



4A I A ) \ A 

where j cu is an external current 



3.1 Electric solutions 



Let us now investigate the Coulomb problem i.e. we will find field configu- 
ration generated by an external, static, point-like electric source: 



J 



Moreover, for simplicity, the source is set to be Abelian. We would like to 
stress that restriction to the Abelian case is not essential. One can easily ana- 
lyze the more general source with three non-zero color components. However, 
the results will be modified only by a multiplicative color-dependent constant 
(cf. jS]). The dependence on spatial coordinates is identical. Because of that 
we remain in the Abelian sector. In order to work only with Abelian degrees 
of freedom one can put 

A°, 



and field equations can be rewritten as follows: 




and 




(9) 



(10) 
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(11) 

Here E ct = —F c0i and E c = E5 c3 i.e. the electric field points in the same 
color direction as the source. We have also assumed spherical symmetry of 
the problem E = E(r)r. The prime stands for differentiation with respect 
to r. Using the Gauss law (JSJ) one can express the electric field in terms of 
the scalar fields 

• ' — " ' I \ ~ b 




A 



(12) 



One can treat this field as the standard Coulomb field in a very non-standard 
medium. Indeed, it can be written as E = ^j- where g e // is an effective 
charge strongly dependent on the scalar fields. Now, we insert it into the 
equations for the scalars and get 

-o-l / i \ -b 

(13) 




and 




-6-1 

(14) 



where a new variable x — 1/r has been introduced. Here dot denotes differ- 
entiation with respect to x. We look for solutions of the equations (|T3)) and 
(|14jl in the power-like form. Namely: 



and 



o(.r) = A\(| + ^) (10, 



1 j,(x)=BA[^ + ^) m , (16) 



,A A 

where constants m, n, A, B and xq are yet to be determined. After some 
easy algebra one can find that 



and 



n = m = (17) 

2+a+b K ' 



{ 4 (a + 6) J • (18) 
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Finally, we have found the following family of solutions of the Coulomb prob- 
lem labeled by the positive parameter (3q 



2 + a + b 



2 + a + b 



The electric field takes the form: 

q 



E 



-A~ a B 



a r>-b 



rA 



2+a + b 



(20) 
(21) 

(22) 



Here we have put Pq = — . Then corresponding energy density reads 

Xq 



1 
2 



e = ~ I q 2 A- a B- b + 



A{A 2 + B 2 ) \ 1_ / J_ 1_ N 
(2 + a + 6) 2 J ^ I WV + J 0/ 



2a + 2b 
' 2 + a + b 



(23) 



Integrating (|2~3j) one finds that the total energy of the fields generated by 
the external static electric charge is finite if a + b > 2 and it is given by the 
formula: 



a + b - 2 2 I y (2 + a + 6) 2 



/5o 



a+6-2 
a + 6+2 



(24) 



Of course, there are also solutions with negative parameter (3q. However, 
such field configurations are not regular for all r. There is a singularity at 
r — > (3. Due to that such solutions are usually interpreted as electric black 
hole solutions with r > j3 jHj- 

In case of the family of solutions with finite energy we can define some 
additional numbers, so-called dilaton (scalar) and gauge ('electric') charge, 
which characterize asymptotic behavior of the solutions. Using the standard 
definitions we obtain: 



Q = r 2 E(r) 

for the electric charge and: 

n 2 d( P 
D * = - r Tr 

dr 



2a + 2b 

qA~ a B- b (3£ +a+b , 



a + b 



D 



2 + a + b 
2 



A(3n +a+b 



2 + a + b 



B(3 t 



a + b 
2 + a + b 



(25) 

(26) 
(27) 



for the scalar charges. 

One can notice that scalar charges and gauge charge are not independent. 
They are connected by the simple relation 

D+D+ _ g a + b_ 



Q Vab 

The finite energy solutions of the Coulomb problem can be interpreted as 
screening configurations. The fields generated by fixed electric charge can 
have arbitrarily small energy. This phenomenon is known from the standard 
non-Abelian Yang-Mills theory where non-Abelian contents of the gauge 
field lowers the total energy Here it occurs even in the Abelian part 
of the model. This suggests that scalar fields can probably represent the 
non-Abelian part of the gluonic sector of QCD [Til] . 

In addition to the family of the finite energy field configurations, there is 
a unique infinite energy solution 

, 2 
/ 1 \ 2+a+b 

0(r) = AA (— J , (29) 

2 

/IN 2+a+b 

ip{r) — BA ( — 7- J (30) 



and 



rA 



, , 1 \ 2+a+b 

E = qA- a B- b A* - . (31) 



rA 



Then the electric potential reads 



(7, -\~ k ~\- 2 , / 1 \ 2 + a + b 

H a + b-2 \rA) v ; 

for a + b ^ 2 and 

U = qAA~ a B- b \nrA (33) 
for a + b = 2. The pertinent energy density is given as 



(34) 

Obviously, the corresponding total energy is infinite. However, for a + b > 2 
the infiniteness of the total energy has its origin in the long range behavior 
of the fields. It is unlikely the standard Coulomb potential in the classical 
electrodynamics where the energy is infinite because of the singularity of 
the potential at r = 0. This effect has been used in many QCD motivated 




7 



theories to model confinement of quarks at the classical level. The confining 
electric potential obtained here does not diverge stronger than linearly with 
distance r. That is in agreement with the famous Seiler constrains for quark- 
anti-quark potential One can observe that the standard linear potential 
emerges only in the limit a + b — > oo which cannot be implemented at the 
Lagrangian level. This means that the model with two scalars does not 
describe the linear potential. However, it can be approximated with arbitrary 
accuracy by taking sufficiently large value of the parameter a and b. In spite 
of the fact that our model does not possess the linear potential, it can still 
be physically interesting. In fact, it gives confining potentials which can 
be compared with many phenomenological potentials obtained from fits to 
charmonium and bottomium states. For example, for a + b = 6 one obtains 
confining part of the Zalewski-Motyka potential 

U ZM = C 1 (yt-^j, (35) 

1 3 

where G\ ~ 0.71Gev2 and C2 — 0.46Gev2. Analogously, the Martin poten- 
tial [T3] is reproduced for a + b = 

As long as the asymptotic values of the scalar fields are not fixed i.e. the 
model does not contain a potential term for scalars, the confining and finite 
energy solutions appear simultaneously. In order to get rid of the screening 
field configurations and to preserve the confining solution one has to add a 
particular potential therm. Let us, for instance, choose it in the following 
form, which enables us to obtain analytical solutions 




nw)=A 4 L«i j [l) ' (36) 

where every parameter fa fulfills 0<fa<4: + a + b. Constants aj are 
restricted by the condition that the potential must be positively valued and 
has global minimum at <ft — ip — 0. Then the equations for the scalar fields 
(JEJ), (fUjl can be rewritten in the form 

+ ^X>(4 + a + &-A) I f) (37) 




and 
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The Gauss law remains unchanged and electric field can be expressed by the 
scalars via the formula In the simplest case, when N — 1 the solutions 
read 

, 2 
/ 1 \ 2 + a + b 

AA{—) , (39) 



i/j(r) 



rA 



, 1 \ 2+a + b 



and 



E = qA 2 A- a B- 



rA 



2+a + b 



(40) 



(41) 



where constants A, B are given by the following set of algebraic equations 
2(a + b) 



a + b) 



aq 



A- a - l B- b + a(A + a + b- j3)A 3+a+b ~ p B , 



and 



{2 + a + b) 2 2 H 



(42) 



(43) 



One can notice that taking into account the potential (|36|) does not influence 
the functional dependence of the confining solutions. Modification of the 
Lagrangian is visible only in the constants A and B. Because of the fact that 
the potential (f3T)J) has unique minimum for the vanishing scalar fields the 
finite energy solutions ()20|) - (f2*2*|) (asymptotically non-zero) can no longer have 
finite energy. Moreover, one can check that such solutions are in contradiction 
with the equations of motion. In other words, the finite energy configurations 
disappear from our model. 

3.2 Magnetic monopoles 

In this subsection we find field solutions generated by the static magnetic 
monopole. In order to do it we adopt the well-known SU(2) monopole ansatz 



A" 



x 



e C ik^(g - 1), A c c 



0. 



(44) 



where g = g(r) is an unknown function. Then the non-Abelian field equations 
read 



4>\ /V 



A 



and 



A) 

a 

' 2A 



9 



0. 



2g' 2 (g 2 - l) 2 

n a 



A 



(45) 



(46) 
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b 

2A 



~2</ 2 , (g 2 - If 



A 



b-l 



Equation (J43|) possesses the trivial solution 

9 = 0. 



(47) 



(4£ 



In other words, we have found the Wu-Yang SU(2) monopole. In the classi- 
cal Yang-Mills theory this monopole has infinite energy due to the singular 
behavior of the gauge field in the vicinity of the monopole. As we will see 
the scalar fields are able to 'regularize' monopole solution. 

After substituting (|4*Hj) to the remaining field equations we can obtain the 
following family of solutions 



and 



CA 



•0 = DA 



J_ 1_ 

rA + J 



rA 



A), 



where the constants read 
C = a^^b 

and 



' b \2 { 2-a- b) / g 2( 2 _ a _fr)2\ 2-a- 

a) { 4(a + b) ) 



D 



\ 4(a + 6) ; 



The energy density originating from (}4*%j) -()50 |) : 



1 



- I C a D b + 



(2-a-b) 2 )^\Va + ]3~o / 



2a + 2b 
2-a-b 



(49) 



(50) 



(51) 



(52) 



(53) 



As one could expect, for a + b > 2 the total energy for the field generated by 
the Wu-Yang monopole is finite and reads 



£ 



a + b ~ 21 - C"D» 
a+b+22\ 



4(C 2 + D 2 )' 
(2-a-b)\ 



Pi 



a + b+2 
a + b-2 



(54) 



Similarly to the case of the electric solution discussed above, the scalar fields 
surrounding the magnetic monopole possess dilaton charges, which take the 
values: 

= — rCpf^ 2 (55) 



a — b 
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For negative value of the parameter /3 in the magnetic solutions (JSOJ) 
we find the magnetic counterpart of the electric black hole configuration 
presented previously. 

It is widely accepted in the literature that finite energy magnetic solutions 
correspond to glueballs (see e.g. |14j ) - effective particles in the low energy 
sector of gluodynamics. As we have shown the model provides the whole 
family of glueball-like states. 

In addition there is a solution with infinite energy: 

2 

0(r) = CA (-) ^ (57) 

2 

^(r) = J DA(— ) 2 "°" b (58) 

for a + b ^ 2. One can check that, for a + b > 2 energy density diverges as 
~ r s at the spatial infinity 

Here 5 takes values from (a + b —>■ oo) to infinity (when a + 6 — > 2). 

Interestingly enough the electric as well as magnetic solutions of the model 
P|) revel an universal property. The solutions describing scalars are (up 
to a multiplicative constant) identical, even though these fields differently 
couple to the gauge field, i.e. with a ^ b. The energy density of the fields 
generated by the electric (magnetic) source behaves as if there has been only 
one (effective) scalar field a coupled to F c ^ u F^ y invariant by a a+b . 

To summarize, there are three sectors of solutions for the model defined 
by (jlj): screening, confining and glueball-like. They appear simultaneously 
i.e. for fixed parameter of the model one can derive finite as well as infinite 
energy solution. Because of the fact that the main aim of such models is to 
provide effective description of the low energy features of QCD it is important 
to get rid of these non-physical screening solutions. It can be easily done by 
adding of the potential term (}36|) . Unfortunately, it removes not only the 
finite energy electric solutions but also all the glueball-like solutions. This is 
highly unsatisfactory. The good candidate for an effective model describing 
the low energy QCD should model confinement and possess, at least one, 
glueball solution. 
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4 The massive scalar fields 



In this section we turn to massive scalar fields. In effective models for the low 
energy gluodynamics such massive scalars are usually interpreted as scalar 
(effective) glueball or/and meson fields. It is unlikely in the massless case 
where scalar fields do not have any clear particle interpretation. Then the 
Lagrangian ([TJ takes the following form 




A I V + 



(60) 

where m^, are masses of the scalars. Let us now consider the Coulomb 
problem for the massive model (JBTlj) and compare it with massless solutions. 
The pertinent equations of motion can be written as follows 




4:7iqS(r) 



and 



V r 2 ^ 




mlil). 



(61) 



(62) 



(63) 



Obviously, the additional massive term for scalars does not inflect the Gauss 
law. Thus, the electric field is given in terms of the scalar fields in the 
following way 

"or 




(64) 



Then remaining equations read 



aq 



-a-l 



-b 



2r 4 A 



(65) 



and 



bq 2 



2r 4 A \ A 



-6-1 



(66) 



Unfortunately, we cannot solve this set of differential equations analytically. 
However, because of the fact that model (jlj) is usually considered in the 
connection with low energy QCD we are mainly interested in the long range 
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behavior of the solutions. In fact, the asymptotic form of the solution for 
r — ► oo is found to be 

_ 4 

<Kr) = GJ-) (67) 



and 



where constants read 



if)(r) = G^ 



1 \ 2+a + b 



(6* 



G A = A 



2m2 yam^ 



(69) 



and 



G,, 



+b- 



-i h 2 



a mt 



2m\ \bm\ 



(70) 



Using we easily obtain the resulting electric field 



E(r) 



/~i—aft—b 

' A- a ~ b 



2+a + b 



(71) 



and the corresponding electric potential 



U(r) = q 



2 + a + b 



(~i—a/~i—h 

^4, 



3(a + b) - 2 A- 



2 + a + 6 



(72) 



One can repeat the previous calculations and conclude that this configuration 
of fields has infinite energy. For a + b > | the total energy diverges due to the 
behavior of the fields at the spatial infinity. Thus, the model (J1J simulates the 
confinement of external electric source, where the confining potential is given 
by f72|). As it was mentioned before, this potential should not depend on 
r stronger than linearly. This provides an upper bound for our parameters: 
a + b < 2. Finally, the massive scalar model describes confinement for the 
following parameters 



a + b G 



3' 



(73) 



In the very special case, when a + b = 2 and masses of the scalars are identical 
rr& = mL = m 2 , we are able to find the generalized Dick analytical solution 



rvv2m \ V 2 m, 




(74) 
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and 

r \a J y v 2m V V 2 m/ 

In agreement with (|T3|) such form of the scalar fields guarantees the linear 
dependence of the electric potential at the large distance. In fact, the electric 
potential reads 

^&)M e2 "- 1+ ^# (76) 

and for r — > oo it diverges linearly. 

In the paper jHj Dick and Fulchert have proposed a model where the single 
scalar field representing the lightest ++ glueball is linearly coupled to the 
gauge fields i.e. <pF* F aiW . It corresponds to a = 1 and b = in our case. 
This particular form of the coupling has been chosen in analogy to the chiral 
quark model ^3]. The confining potential derived from their model reads 
ra. However, one can notice there is another particle which is relevant in 
the low energy regime - a scalar meson. It seems to be reasonable to assume 
that both scalars inflect the low energy dynamics. Due to that one should 
consider these scalar fields together. Then the natural generalization of the 
Dick-Fulchert model is a theory with one glueball-meson coupling, that is 
<j)ijjF® u F aflu (a=l, b=l in our model). It is very interesting that this form 
of interaction ensures also the linear electric potential. This unexpected and 
striking result can suggest that the correct low energy effective action should 
consist of more than only one scalar field. 



5 The non- dynamical scalar fields 

In the previous sections we observed that the dynamical scalar fields (massive 
as well as massless) can strongly modify electric and magnetic solutions. 
Now, we will show that these fields can play very non-trivial role even when 
we would drop the kinetic term for scalar [16J. As it will be clarified later a 
model with non-dynamical scalar fields can be also interesting in the context 
of the low energy gluodynamics. 

Let us start with the mixed case: dynamical field <fi and non-dynamical 
ip. Additionally, we assume that the potential for the non-dynamical field is 
analogous to that in the formula (jHfij) . Then the Lagrangian (in the simplest 
version) reads 
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where d is a positive parameter. We will analyze the Coulomb problem. 
Then the pertinent equations of motion take the following form 



E 



4nq5(r) 



(78) 



and 



V 



3 / V> 



aE 2 U 

d—l , 17,2 



adA 6 I j 



bE 2 uy f 

2A UJ U 



6-1 



(79) 



10) 



Obviously, the variation with respect to the non-dynamical field gives a con- 
strain on the scalar field ip (J80J) . We can solve it and express ip in terms of 
the dynamical fields (f) and E 



't\ = (- 



bE 2 



i 

d-b 



v A J \2adA 4 \A / 
Then the remaining field equations 



r 2 I — 



ad 
d-b 



b -. / 

d-b 



( bE 2 \ 
A J \2adA 4 ) 



(82) 



and 



aE< 



f bE 2 \ 



--1 



2A \2adA i ) I A 



Finally, we obtain the following solutions 



id 

I \ ad+2(d+b) 

(/) = /-/A(4 

.rA, 



and 



E(r) = A 2 qd+b 



. L. 4(d-b) 

| // — f 2 " 1 " : 

2ad/ VrA, 



Here constant if takes the form 



b+a 



2d ( b \ / a (ad + 2(d + b)) 2 \ - d +^ d + b ) 

JJ = gad+2(d+b) III I 



2ad 



A(d-b)(ad + Ab)^ 



(83) 



34) 



55) 



(86) 
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The corresponding electric potential reads 



ad-2(d+b) 

<t± ad + 2(d + b) I b \ d+b ^-juL ( 1 \^+2(d+b) 



It has confining-like behavior if the parameters fulfill the following condition 

l> 2( ^ + n"1 >0- (88) 
2{d + b) + ad K > 

On the other hand, because of the fact that the non-dynamical field play 
role of the Lagrange multiplier, the action can be rewritten in terms of the 
dynamical fields only. The new Lagrangian has the form 

ad fj 

1 ( d>\ d - b /F\~b _____ i 



L = -i[i) [a) FC ^ F ^ + 2 d ^ d ^ ( 89 ) 

where F = \F C » V F* V . One can easily check that, in the case of Coulomb 
problem, the pertinent equations of motion are identical to (J%2|) and (J%3|) . 
The system with the dielectric function governed by dynamical and non- 
dynamical scalar field is equivalent to the model where the permittivity de- 
pends on the dynamical scalar as well as gauge fields. Thus, we have obtained 
a generalized version of an effective theory proposed by Pagels and Tomboulis 
a long time ago |17j . On the contrary to the original version (|%9*|) contains 
non-minimal gauge-scalar coupling along with an exponent of the standard 
Yang-Mills invariant. This equivalency seems to be particularly interest- 
ing since the Pagels-Tomboulis model (and the similar Adler-Savvidy model 
|18j . |19j . pU]). reproduces the trace anomaly already on the classical level 
and gives correct prediction for the confining potential. 

Obviously, it is possible to obtain the standard Pagels-Tomboulis model 
treating also the second scalar field as non-dynamical. Let us assume that 
the potential for this field is 



V(4>) = a' 



A 



where e is a new parameter and a' is a dimensionless constant. Then both 
scalars can be written as functions of the electric field. Finally, we obtain 

where 

ad + be 

2o 



ed — ad — be 
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As it was shown in J7|, [21], this model guarantees the confinement of ex- 
ternal, Abelian electric sources for 

5>\. (91) 

Duality demonstrated above allows us to believe that an ultimate gluodynam- 
ics effective theory would share properties correctly described by all instances 
presented in this paper. 

6 Summary 

The main aim of this work has been analyzing of electric and magnetic so- 
lutions in the standard classical Yang-Mills theory coupled to two scalar 
fields. For particular values of the parameters of the discussed models (with 
massless, massive as well as non-dynamical scalar fields) we have obtained 
confining-like field configurations. External, static, electric sources have in- 
finite energy due to the long range behavior of the fields. From this point of 
view all models are very similar. Clearly, the values of the parameters are 
model dependent and different potentials known from fits to the phenomeno- 
logical data are derived for different a and h. One can also notice that the 
model with massless scalars is not able to describe the linear electric potential 
which corresponds to the limit a,b — > oo. In the other words this model can 
only approximate the linear potential (but with arbitrary accuracy). More- 
over, in contradiction to other Lagrangians, there are finite energy electric 
solutions in the massless scalar model. Here, such solutions representing 
screening phenomena at the classical level, are rather non-physical. They 
are removed from the spectrum of the theory by adding a potential term for 
the massless fields. 

In spite of the fact that presented models in the similar way describe the 
confining solutions they differ profoundly. 

The main difference between our models has its source in the physical 
interpretation of the scalars. This is strongly connected with the glueball 
problem. Let us firstly discuss the simplest case i.e. model with the massive 
scalar fields. Then these field can be identified with some massive parti- 
cles relevant to the low energy sector of QCD. Evident candidates are the 
scalar glueball ++ and a scalar meson. Both objects exist due to the non- 
perturbative effect, confinement of the gauge and quark fields respectively. 
Following that one could expect that these field should appear in the 'demo- 
cratic' manner. They should couple with the gauge field identically: a = b. 
This gives us prescription how correctly generalize QCD motivated Dick- 
Fulcher model. Namely, instead of one glueball coupling we introduce one 
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glueball-meson coupling. Surprisingly, such form of the dielectric function 
provides the standard linear electric potential. For identical masses of the 
scalar glueball and meson the analytical solution has been explicitly found. 

In case of the massless fields situation is a little bit more subtle. Mass- 
less scalars cannot be rather identified with any particles. There are not 
massless particles in the low energy sector of QCD. Even thought scalars do 
not possess particle-like interpretation the existence of the glueballs is not 
excluded. There are at least two possible ways of introducing glueballs in the 
framework of massless scalar model. Firstly, they can appear as finite energy, 
magnetic solutions. However, such configurations exist simultaneously with 
electric screening solutions. If we get rid of the screening solutions by adding 
a potential term to the Lagrangian then also finite energy, magnetic solutions 
disappear. Unfortunately, we are not able to preserve this glueball-like sector 
and remove the screening sector. On the other hand one could try to find 
glueballs spectrum in the analogous way as it has been done for charmonium 
and bottomium states i.e. by means of Schrodinger equation. In this pic- 
ture glueball would consist of magnetic monopole and anti-monopole. Then 
the potential in Schrodinger equation would be potential between magnetic 
monopoles |22j- Similar glueball model has been recently considered in the 
paper [2H]. 

Even more interesting from glueball point of view is the model with the 
non-dynamical scalar field (or fields). Such model corresponds to the gen- 
eralized Pagels-Tomboulis theory. Here glueballs could appear as toroidal 
solitons with nontrivial Hopf index [23]. It follows from the observation that 
the restricted version of Pagels-Tomboulis model with 5 = — \ possesses such 
solutions |2Sj (but then electric sources are not confined). We expect that 
additional scalar field could joint the glueball and confining sectors. 

To conclude, all models presented here can be treated as candidates for 
the effective model for the low energy QCD. They describe confinement and 
give reasonable quark-anti-quark potentials. However, the problem of the 
glueballs has not been satisfactorily solved yet. We would like to analyze it 
in our forthcoming paper. 
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